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ABSTRACT 


This thesis is a design of an alternative system which may provide heating to the Naval 
Postgraduate School swimming pool. Particularly, it is a solar heating / photovoltaic 
system designed for a better efficiency and less cost of installation and maintenance. 
Principles of heat transfer, control and fluid dynamics theory are used for the determination 


of this heating system elements. The feasibility of its installation and use is analyzed. 
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I. INTRODUCTION 


A. GENERAL 

Nowadays, high energy consumption based on fuel has been associated with 
economic development and a high quality of life. With time, energy users’ demand 
increases, while the fuel reserves are dwindling. The conflict between demand and supply 
causes an escalation of the fuel cost, which should mandate the research, development and 
use of alternative resources of energy. We should develop alternative means of energy and 


fuel conservation. 


B. OBJECTIVES 
As an example of a heating usage needing optimization, we will analyze and design a 


new heating system for the swimming pool of the Naval Postgraduate School. 


C. SYSTEM DESCRIPTION 

The existing system, using fuel to heat the water and electric power to drive the pump 
for circulation through the heat exchanger and the filters, is old, wastes fuel and is 
inefficient. Presently, the boiler plant is used to keep the pool at the desired temperature. 
New solar system designs here offer several alternative energy schemes with various 
tradeoffs. An optimum design is suggested. 

The objective is the design and modeling of a system using a minimum area of 
collectors and electric power to provide a maximum thermal energy to the pool. Taking 
into account its geometry and masonry materials, as well as accurate weather data of the 
Monterey area taken from meteorology statistics (Table 1), and applying heat transfer 
principles, a total heat loss is calculated. That heat loss is balanced by the collectors and the 


solar irradiation providing thermal energy to the water and maintaining its temperature. The 











heat loss due to conduction, convection, radiation and evaporation is calculated for each 
month in Chapter II. 
Power vs cost tradeoffs are made for low-power DC pumps with photovoltaic and a 


separate filter pump system and standard AC. 


TABLE 1. WEATHER DATA FOR THE MONTEREY AREA. 







RELATIVE 





MONTH MONTHLY TEMPERA- _ (degrees F) 





HUMIDITY 





TURE 


(mph) “a 

maximum minimum 

june | oro | sor | ses | om | 
uy | sso | sis | 97 | scm | 9 | 
august | 0 | 526 | coe | ses | 79 
foctoner_| 703 | soo | coe | seis | 76 
peemee | st_{_ anf ssf ccs | 
[pecember | oro | x7 | sas | sci | os 





































II. HEAT TRANSFER ANALYSIS 


A. HEAT TRANSFER DUE TO CONDUCTION 

It is considered that the heat is transferred by conduction under one-dimensional, 
steady-state conditions. The one-dimensional refers to the fact that one coordinate is 
needed to describe the spatial variation of the temperature, and the heat transfer occurs in 
that direction only, as can be seen in Figure 1. The system (ground-walls of the swimming 
pool water) is characterized by steady conditions, since there is no energy generated per 
unit volume of the medium of the wall, assuming that the temperature is independent of 


time. The heat diffusion equation is, 


where, 
q=energy rate, which is generated per unit volume of the medium 


os : 
OS a, = time rate of change of the internal (thermal) energy of the medium 
t 
per unit volume. 


can be expressed as: 


The general solution of the eq.(2.2) is: T(x) =C,x + C,, with T(0) = T, , and T(L) = 
Tg». Applying these conditions to the equation we have : C,= Ty, , and atx =L, Ty 5= 


C,bL+C,=C,L+ Ts 1 that means C, = (Ts 4 : Ts, )/L. 








Figure 1. Heat transfer through a plane wall. 


Substituting into the general solution, the temperature distribution is then : 


TO? (Tes Teale tle, 


It is evident from this result that the temperature varies linearly with x and its rate of 


change with respect to the x coordinate is: 
dT _ 1 
ox 7("s,27' si) T (2.3) 


Then the conduction heat transfer rate is obtained by using the Fourier’s law equation, 


dT kA kA 
Ge eae a 8 SA ea Se 











In our case having a series composite wall (Figure 2), the temperature distribution is 


the sum of the distributions for each meu.um separately. Therefore, the temperature 


distributions for the insulating fiber resin and the concrete block walls are T(x) and T.(x) 


respectively. 
xX 
a ca cad a a (2.5) 
and 
ww = oe 
Tax =Cl T) i ok, (2.6) 


where, T, is the temperature at the thermal contact between the two walls. 

Taking the thermal conductivity coefficients of the fiber resin (k,= 0.038 W/mK) and 
the concrete block (k,= 0.67 W/mK) from the Tables A-1 and A-2 respectively, the 
thicknesses of the walls (L=0.01m and L,=1.524m), the temperatures of the water, and 
the ground (T, and T,) and applying into the eq. (2.4), the inlet and outlet transfer rates are 


obtained. 





Ground 





Figure 2. Temperature distribution into the swimming pool walls 


The inlet heat rate q, ,,, which is transferred through the composite wall, from the hot 
water (swimming pool water at 80°F=300°K) towards the cold ground, which is assumed 


to have a temperature 2°K lower than the mean temperature of the air, is equal to the outlet 


heat rate q, out 


k,A kA 
Gy, in =, on L, or aes oy Pai a ad 
k(T,- T k(T, -—T 
rs Me 7 u ott 7 e) aes 


f 





Plugging in the values in the above equation we obtain the temperature T, and then 
from the eq. (2.7) or (2.8), the conduction heat transfer rate q, 18 calculated for each wall 
and the bottom. Summing up the heat transfer rates for each particular area A,. the total 
heat loss due to conduction is obtained for January (Table 2). Similarly working for the 
mean temperatures in Table 1, the conduction heat losses for each month are calculated, as 


shown in Table 3. 


TABLE 2. CONDUCTION HEAT LOSS (JANUARY). 


F_fayeasn’s |g osm) [agccm’) | ag cam Atom 
331.05 | 478.19 | 161.85 | 478.19 | 4796.64 






Therefore, q,,,g= 6245.93 Watt = 15.86 MBtu per month (for January). 





TABLE 3. MONTHLY CONDUCTION HEAT LOSSES. 


MONTH MB1u per month 
£245.93 
5780.9 

l 


3.26 





















: 3 
5205.5 


5931.46 


B. HEAT TRANSFER DUE TO CONVECTION 

The term convection refers to energy transfer which occurs between a fluid in motion and a 
bounding surface when these two are in different temperatures. In this study the moving 
fluid is considered to be the atmospheric wind blowing in average speeds as in Table 1 for 
each month, and the bounding surface is the surface of the water of the pool. Since we 


have a flow caused by external means, we speak of forced convection. 





Considering a small area A, of the water surface with temperature T,, the wind speed 


V, and air temperature T..; a local heat flux q occurs which is proportional to the 


difference (T, — T..), that is: 
q’=h(T,-T..) 


where, h=local convection coefficient. 


By integrating the local heat flux q over the area A,, the total heat transfer rate q is 


obtained. 


q=), 4A, 


and 


q=(T,-T.) fy hdA, 


Then defining the average convection coefficient h for the entire surface, the total heat 


transfer rate is expressed as: 


q=hA,(T,-T..) (2.10) 


Setting the water surface area A,= A = 652 m?, the water temperature 5 OF =T,, = 80°F = 


300°K and the air temperature T = T,= 56°F = 286.33° K (for May) the eq. (2.10) 


becomes: 


q=h A(T, -T,) (2.11) 














Figure 3. Heat transfer from the surface of the water to the air due to the 
forced convection. 


Assuming that the wind is blowing lengthwise with an average speed as shown in 


Table 1, we calculate the Reynold's number, 





where, V is the wind speed, L the length of the swimming pool, and v,., the kinematic 
viscosity of the air. From the Table A-3, for the mean air temperature T, of each month we 


get v,., the conductivity coefficient k of the air and the Prandt! number Pr. Thereafter, 


air’ 


calculating the Reynold's numbers for each case, we see that the flow over the water 


surface is turbulent, because 5x10°<R,<108. Since 0.6<Pr<60 the equation 


Nu=(0.037R,49-871)Pr'8 (2.12) 








is used (Ref. 3] for the calculation of the Nusselt number (Nu), which provides a measure 
of the convection heat transfer occurring at the surface. The Reynold and Prandtl numbers 
are applied to eq. (2.12) and we get Nu. As described in [Ref. 3:pp.337], the defining 


expression of Nusselt number is: 


hL 


Nu= ars (2.13) 


Solving eq. (2.13) with respect to h' and substituting the obtained values of k, L, and 


Nu, we get the convection coefficient h’. Finally, the amount of heat transfer due to 


convection is calculated by applying the values of h’, A,T,, and T, to the eq. (2.11). Then, 


we have the results for each month as shown in Table 4. 
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TABLE 4. MONTHLY CONVECTION HEAT LOSSES. 














C. HEAT LOSS DUE TO THERMAL RADIATION 

The heat transfer by thermal radiation is an important process which occurs at a single 
surface, in this situation, and depends strongly on the surface geometry and orientation, as 
well as on the water radiative properties and temperature. Thermal radiation is related to 
energy released as a result of oscillations or transitions of the electrons which constitute the 
water. These oscillations are sustained by the internal energy, and therefore, the 


temperature of this matter. The water of the pool, being in a higher temperature qT, than 


that of the surrounding air, tends to achieve a thermal equilibrium with the atmosphere by 


i 








getting cooler. This cooling is associated with a reduction of the internal energy stored by 
it, and the consequence is the emission of thermal radiation from the surface. This 
radiation is the energy rate emitted by the water as a result of its temperature T, and goes on 
until the matter reaches the atmosphere temperature. 

The nature of the thermal radiation is a propagation of electromagnetic waves and more 
particularly those with a wavelength from approximately 0.1 to 100 um. This portion of 
electromagnetic radiation spectrum includes a part of ultra-violet, all of the visible and the 
infrared. The rate at which the radiation leaves the surface of the water is determined by the 


equation: 


dus=AFoe(T, - Tan ) (2.14) 


as in [Ref. 3, pp. 760], where, 

A = Area of the surface 

F = View factor 

© = Stefan-Boltzmann constant 

€ = Emissivity of the water = 0.96 (from Table A-5) 

T, = Water temperature 

Tum= Atmosphere temperature = 273°K. 

The view factor is a dimensionless number which is defined as the fraction of the 
radiation leaving the surface of the water and is intercepted by the atmosphere. For this 


study the view factor is considered F = 1. However, for each month the amount of heat 


transfer by thermal radiation is the same. 





Grad = 652m? x 1x 5.67 x 10° W / m?- K* x 0.96 x (300* - 273*) = 90336 Watt 


Then the monthly heat loss due to radiation is: q,, = 221.9 MBtu per month. 


D. HEAT TRANSFER DUE TO EVAPORATION 

The evaporation from a free-water surface is a process associated with water-vapor 
transfer from the water surface to the air and occurs by the combined processes of 
convection and diffusion, Figure 4 is the schematic model. 

The surface of the pool at temperature T, is exposed to a stream of moist air. Into the 
boundary layer, the velocity of the air V(y) varies from zero to V, at the outer edge. 
Temperature T(y) decreases from T, to T,, and the humidity decreases in the y direction 


from %, at the water surface to >, in the air. The transfer of mass of water-vapor having an 


internal energy causes a thermal loss whose rate is given by the equation: 


evap = Amb (1 — 02)P a, aA (2.15) 


where, 
h_, = mass transfer coefficient (m/s) 
he = heat of vaporization 
o, = relative humidity 


P, sar = density of the saturated air 


A = area of the swimming pool. 





Water surface 





Figure 4. Evaporation of water into moist air. 


The mass transfer coefficient is given by the equation: 


, 


h 


hy =——_— 
p.c, (Sc /Pr) 


where, 
h = convection coefficient 
Pp, = air density 


Cat specific heat of the air 


Sc = Schmidt number 


Pr = Prandtl number. 


(2.16) 





The values of the above variables in eq. (2.15) and (2.16) are taken from [Ref. 3 
Tables A-3 and A-6} and correspond to the mean temperatures of each month. The relative 


humidity of the air ®, is the mean value for each month taken from Table 1. 


However, the heat transfer due to evaporation is calculated and shown in Table 5. 


TABLE 5. MONTHLY EVAPORATION HEAT LOSSES. 
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Ill. SYSTEM OPERATION 


A. SOLAR INSOLATION 

The average solar radiation density outside the atmosphere is approximately 1395 
W/m? [Ref. 6] for a plane normal to the rays of the sun. This solar constant, reaching the 
earih’s surface, reduces because of the atmospheric reflectivity and absorption, which 
depend on the angle of the incident radiation as well as the local climatic conditions. In 
general, the average reflectivity is 35% and the absorption 20%. Then the average incident 
power density is approximately S = (1-0.35)(1-0.2) 1395 W/m? = 725 W/m? and it occurs 
at noontime (average maximum value). Designing a solar heating system is of interest to 
calculate the daily and then the monthly average incident power density, which varies with 
the time (seasonal variation) as well as the place (geographic position). 

1. Geometry of the Incident Solar Radiation 

The incident solar power varies due to the daily rotation of the earth and to the 

varying tilt of its axis with respect to the equatorial plane. This angle varies for a given 
place in a sinusoidal sense (Figure 5) through the year and gets its maximum value on June 


22 (23.59) for the northern hemisphere and its minimum one on December 22 (-23.5°). 


The above angle is called declination of the sun ®,. Considering a spherical coordinate 
system (Figure 6) with its origin to the center of the earth, where a, is the normal vector to 
the surface of a specific place M, 1/2-6 is its latitude and $+n/2 is its hour angle, the 


incident rays of a flux S impinge on this surface with an angle 8, with respect to the 


equatorial plane. In this case, for convenience, it is assumed the sun rays to be in the y - z 


plane, so noon coincides with 0 = 77/2. 


17 








Declination (degrees) 

















Figure 5. Declination of the sun through the year. 


Then the normal incident flux is: 


S,=-a,S (3.1) 


Resolving the vectors a, and S into rectangular components, we have: 


a, = a, sin@ cosd + a, sin Osing + a, cos@ 


S= “ay S cos8,; - a, S sin, 


and substituting in eq. (3.1), we get: 


S, = S cos, sin8 sing + S sin®, cos® (3.2) 








Figure 6. Polar coordinate system for solar radiation. 


The above equation gives the incident solar radiation with respect to time (@ and 


8.) and position (7/2-8). To obtain the average daily power density for a particular day we 


integrate it over the daylight time, that is: 


J (cos 6, sin @ sin > + sin @, cos 8)do (3.3) 


daylight 


The corresponding hour angles for the sunrise and sunset are 9, and (1-9,,), 


respectively. The sunrise angle 6,, corresponds to that 6 which makes the integrand of eq. 


(3.3) zero, that is: 








cos@, sinO@sing,, + sinO, cos® = 0 








: _ tan 04 
Om = tan@ 
rr tan 0, 
ds a tan 8 On) 
And the eq. (3.3) becomes: 
Seco ~f (cosOq sinOsind + sin@g cosO)do 
2 
= a (cos64 sin Osing + sin@y cos@)do 
18 a 
= = (00564 sin@cosom + sin Oy cos0(x/2-m)| (3.5) 


In the above eq. (3.4), using the mentioned solar constant (725 W/m?), and 
substituting 6 = 53° 25 (Monterey latitude = 36° 35’ N), as @, the mean declination of the 
sun for each month taken from [Ref.10] and as daylight period the mean for each month, 
the average daily solar insolation is obtained, as shown in Table 6. These results are 
consistent to the data based on weather statistics and given in [Ref. 7]. 

Zz. Optimization of the Solar Collector Tilting 

In designing a solar heating system, the solar collector tilting is taken under 

consideration in order to maximize the received solar power density. The basic concept for 


a fixed collector is that it must receive solar energy for most of the daylight. 
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TABLE 6. MEAN DAILY SOLAR INSOLATION (MONTEREY). 


P| angtey rgay | pu reeaay | wim | 
April 
. Po | naz 













March 


September 250 





Assuming that a, is normal to the collector surface and © is the angle between 
the normal and the vertical a,, then the incident radiation on the collector is: 
S,= -a,8 (3.6) 


n 


From the geometric relations (Figure 7) substituting in eq. (3.2) the angle 8 by 


(8 + @), we get: 
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Figure 7. Tilted collector at an angle a with respect to the ground. 


S,, = 5 cos, sin(6+a) sind + S sinO, cos(0+x) (3.7) 


Integrating eq. (3.7) over the daylight period, the average solar irradiation S_, ,, 


is obtained. 


J [cos 0, sin(@+a)sin > + sin 8, cos(O+ a) ]do 


daylight 


[ [cos 8, sin(6 + asin 9 +sin 6, cos(O + a)] do 
§ ,t/2 . 
= + \. [cos 8, sin(@+ a)sin @ +sin 8, cos(O + a)]do 


— feos 6, sin (8 + a)cos >,, + sin 6, cos(@ + a) (> - om) | (3.8) 
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In order to find the optimal tilting & of the collector for the winter (negative 
values of sun's declination), we define a variable S_ ,,,,; which is the sum of the average 


solar irradiation S_,, for each day over the winter time (Sep. 22 through Mar. 22). 


Then S n total = we S n av(% 84.9 m) 
int er 


The total solar irradiation is a function of the angle a and a specific value of it 


maximizes S Numerically, we obtain the angle @ after a series of runs over a range 


n total’ 


of feasible angles and record the resulting S_, ,o.4;. The change of the S_ |.) during winter 


with respect to the tilting angle a, is shown in Figure 8. It can be seen that S_, ,,.,) 1S 


maximized for a = 56°. 
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Figure 8. Total solar irradiation received in winter (Sep. 22 through Mar. 
22) vs tilting angle of the collector. 
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Figure 9. Total solar irradiation received in summer (Mar. 22 through Sep. 
22) vs tilting angle of the collector. 


If the collector is tilted to operate during the summer (Mar. 22 through Sep.22) 
and have the maximum efficincy, working similarly we obtain the optimal angle a = 10° 


(Figure 9). For a year-round operation, the optimal angle is a = 33° (Figure 10). 
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Figure 10. Total solar irradiation received by the collector through the year 
vs its tilting angle. 


B. ENERGY BALANCE 


In Chapter II the net heat rate loss Q, due to heat transfer is calculated for each month 


(Table 7), taking under consideration the weather conditions, the geometry and the 


construction of the swimming pool. Solar irradiation qg has also been calculated (Table 6) 


for each month and for the Monterey area. 


Solar power Q, is collected by the array and through a heat exchanger delivered to the 
thermal load (Figure 11). Considering the average efficiency (fraction of the solar 


irradiation extracted as useful energy) of a collector usually used in this type of heating 


systems n = 0.7, then the power drawn by the collector is Q. =n Q,. Assuming also that 


the heat exchanger provides to the pool all the thermal energy it receives from the array, and 


additionally its free surface Ap receives directly a part of the solar power. we model the 


collector - pool system. 
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(Thermal load) 


Figure 11. Block diagram of the energy balance in the collector - pool 
system. 


Qe + (1-Ps) ds Ap = Q 
or 


n Qs, + (1-Ps) dg Ap = Q. (3.9) 


where Qg, Qc, Q; are energy rates (Watts), Pg is the water reflectivity taken from [Ref. 3], 


and (1-p,) is the fraction of the solar power received by the swimming pool surface. 


The above relation denotes the energy balance that must exist for the water temperature 


to be maintained. Eq. (3.9) becomes: 


n Acgg + (1-Ps) dg Ap = Q, (3.10) 
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TABLE 7. TOTAL MONTHLY HEAT LOSS. 


MONTH NET HEAT LOSS MONTH NET HEAT LOSS 
(WATTS) (WATTS 
















where dg, q; are power densities (W/m?), Ps = 3 and A¢ is the area of the collector array 


(m7) to be calculated. Designing the heating system to operate all year, we are looking for 
the worst case, where the array will provide energy to balance the greatest loss while the 
incident solar irradiation will be the least (maximum aperture). This turns to be Ac= 2183 
m? for December. If the heating system is to operate from April through October, then Ac 
= 720 m* = 7750 ft’. 

The above calculated areas of the solar collector array, which is to be used in each 
case, are too large and off the swimming pool site. Up to this point, the design of the 
heating system was done without taking under consideration the already existing swimming 
pool cover. Assuming that the pool will be open from 10 a.m. to 6 p.m. and covered the 


rest of the day, then the net heat loss Q, will be: 


Q 7 Qnae + 1/3 Qeanaeeh + 2/3 Qe ond ciy + 1/3 9 + 1/3 Qhad.srf + 2/3 OD aikeny (3.11) 








where, Q.ona.w 1S the conduction heat loss through the walls and the bottom of the pool 


which remains the same after covering (Table 3), Q , 1S the convection heat loss from 


conv,sr 
the free water surface which holds for one third of the day (open pool) (Table 4), Q. ong cov 
is the conduction heat loss through the pool cover which holds for two thirds of the day 


(covered pool), Q,,,,,, 18 the evaporation loss for open pool only (Table 5), Qrig grr iS the 


evap 
radiation loss through the free water surface (one third of the day), and Oey is the 
radiation loss from the pool cover. This cover consists of two polyethylene films and a 
middle cross linked polyethylene foam layer (Figure 12). 

To calculate the conduction heat loss through the cover, we use eq. (2.7) or (2.8) 


taking under consideration as insulating material the polyethyiene foam middle layer only. 


A 
Qconv.cov a —2_(T, 2, T,} 


where ko is the thermal conductivity of the middle layer of the cover, L, is the thickness of 
the layer, ib the temperature of the pool and T,, the average temperature of the air and A the 
area of the pool. Taking the thermal conductivity coefficient of the polyethylene from the 
Table A-3 (k, = 0.015 W/mK) and for L, = 0.005 m, T, = 300° K, T, = 286.33° K (56° 
F for May) and A = 652 m?, we obtain Q 


cond,cov 


= 46944 Watt (for May). 


For the calculation of Q we use eq. (2.14). Then: 


tad,cov’ 


OO aaieav =A Foe(T,* ~ T*am | (3.12) 











Figure 12. Cross section of the swimming pool cover: (a) Polyethylene scrim, 
(b) cross-linked polyethylene foam, (c) polyethylene film. 


where A is the area of the swimming pool surface, F is the view factor and F = 1, o is the 

Stefan-Boltzman constant, € is the emissivity of the outer polyethylene film (€ = 0.85 taken 
= 3000 kK = 27730K; ; 

from Ref. 3), a = 300" K, and T,,_,= 273" K is the atmosphere temperature. Applying 


these values in eq. (3.12) we get Q = 79985.34 Watt. 


rad,cov 
Substituting in eq. (3.11) all the values of the different losses for a particular month 


(e.g., for May) we obtain: 


Qy = 5242.26 + 0.33 x 60754.73 + 0.67 x 46944 + 0.33 x 25615.48 + 


+ 0.33 x 90336 + 0.67 x 79985.34 = 148597.96 Watt. 


Comparing the above result with the total heat loss for May from Table 7 (181948 
Watt), we get energy savings 33350 Watt or 18.32%. Working similarly for all the month 


heat losses, we get an average energy savings of 20%. Then using eq. (3.10), we 


recalculate the area of the collector array A, looking again for the worst case (greatest loss - 
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least solar irradiation), while the pool cover is used to reduce the heat loss of the swimming 
pool during the night hours. For an all-year operation of the heating system, it is required 
an array A, = 1616 m?. If the system operates from April through October, then A, = 446 
m? = 4800 fr. 

For a high performance solar heating application, as in this case, a commonly-used 
collector is the glazed solar collector made of copper plate with copper piping, (Figure 17). 
A better conductivity of the thermal energy is provided by copper fins roll-formed around 
the copper tubes. The collector surface is glazed with a high-quality glass which operates 
as a heat trap, reducing radiative and convective losses emitted by the collector. The 
dimensions of a formal collector panel of this kind are 4'x10'; for the heating system 
operating in the summer, 120 panels are required and their cost is $1300/panel or $156,000 


for the whole collector array including its installation. 


C. PUMP AND PIPING 

The solar heating system includes two loops of fluid circulation. The collector loop 
consists of the solar array, the heat exchanger and a pump P, for circulating the water 
through the collector and back to the exchanger (Figure 13). The second loop circulates the 
water from the exchanger to the filters and then to the pool using a pump P,. This pur) 
already exists in the old heating system. It has a flow rate of 750 gal/min circulating the 
pool water roughly three times a day and is driven by standard ac. The other pump P, is to 
be determined with respect to power (h.p.), that is, we have to determine a flow rate (m*/s) 
for achieving a good conduction of heat from the collector to the water as it passes through. 
On the other hand, a large flow rate through the exchanger is a poor practice resulting in 


excessive pumping and cost. Then P, must be adequate for delivering sufficient water 


flow for good solar collector performance during the summer months where the incident 


solar irradiation is at its highest (300.47 W/m? in June). 








Figure 13. Schematic model of the heating system. 


Assuming that all the heat is transferred from the heat collection fluid to the swimming 


pool water, then the thermal energy in the exchanger is: 


E, =(m,c,)AT, -(m,c, )AT (3.13) 


where m, is the water mass into the collector loop, C, is the specific heat of the water in an 
average temperature 330°K in this loop and c, = 4184 J/kg °K, AT. is the decreasing 
temperature in the exchanger, m, is the water mass into the swimming pool loop. c, is the 


specific heat of the pool water in an average temperature 300°K in this loop and as 4180 
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J/kg °K, and AT, is the increasing temperature in the pool. Dividing both parts of 


eq.(3.13) by time, we get the delivered power by the exchanger to the load. 


Q. = (mh,c, )AT, - (m,c, )AT, (3.14) 


where m, and mp, are the mass flow rates (Kg/s) provided by P, and P, respectively. Eq. 
(3.13) shows that the rate of change of the heat in the exchanger is determined by the rate 
of heat addition from the collector less the rate of heat removal by the swimming pool 


(thermal load). Since all the heat is transferred from the collector to the exchanger we have: 


nA.qs = Q. =(mm,c,)AT, —(m,c, AT, (3.15) 


Determining a control strategy where the collector pump P, is turned on when AT, = 
Tsa= es 4°K or AT a= 2°K below the desired average temperature of 300°K, and 
taking under consideration that M, = 750 gal/min = 47.18 kg/s and q, = 300.47 W/m? (for 
June), we substitute in eq.(3.15) and solving for m, ,we get a mass flow rate 29.17 Kg/s 
or 470 gal/min. 

Assuming that the piping of the heating system has a 200' length of 5" pipe and that 


the fittings (90s, tees, 45s, couplers) have an equivalent length of 100' of pipe , then for a 


flow rate of 500 gal/min, we get from Table A-7 a head loss 3x4.61' = 13.83’. 
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Then: 


Total head loss in pipe and fittings : 13.83' 
Height from ground to roof : 10.00’ 
Approximate friction loss in plumbing: 60.00" 
Total head loss . 83.83' 


With the operating - system characteristics: 
a. Total head H =84ft 
b. Pump discharge U*= 470 gal/min 
we have to determine the pump P, size, revolutions and break horsepower for water 
circulated at an average temperature 330°K. In Figure 21 [Ref. 11] using the measured - 
performance curves for a particular model of a centrifugal water pump with an impeller 
diameter D = 35 in operating in 710 r/min, we get head H = 187.5 ft, discharge U = 19000 
gal/min, and a break horsepower P = 1000 bhp. 

Given the above data for a mode! pump and correlating it with our system 
characteristics we may calculate the break horsepower of P, which will be of the same 
geometric family (centrifugal) and operate at an homologous point (same performance). 
This correlation can be done by using the similarity rules (Ref. 11], that is, head goes as 
n*D2, discharge as nD3, and power as pn>D°, where n is the number of pump revolutions 


and p the water density. The similarity rules are: 


* * 2 . 2 
HL [=| (2) (3.16) 
H n D 


* * «\3 
U -“(2) (3.17) 
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Applying the values H, D, n, taken from the curves, and H", U" for our system, in the 


eqations (3.16), (3.17) we get an impeller diameter for P, D” = 6.73 in and n” = 2472 
r/min. The data of the model pump taken from figure 21 corresponds to a water 
temperature of 60°F or 288°K with density p = 1000 kg/m?. In the collector loop (135°F 
or 330°K) the water average density is p*= 984 kg/m. Substituting p*, n*, D’, P, p, n, 
D in eq. (3.18) we obtain P"= 11 bhp or 8.13 kW. This pump costs approximately 
$3,000. 


D. CONTROL UNIT 

The controls of any solar heating system can be considered its brain. With the aid of 
sensors, it can determine when heat is available for collection and make the decision when 
to start and stop the heat collection process. The control systems used in such solar heating 
systems are called differential temperature controllers. Basically they compare the 
temperature of the solar panels to that of the load, (Figure 14). When the temperature of 
the panel exceeds that of the load by 4°C, then a signal is sent by the unit, and the actuating 
devices are turned on. The collector-to-load differential at which the unit turns on is called 
turn.on delta T. The collector-to-load at which the unit turns off is called turn off delta T. 
The difference between the turn on and turn off delta T is called temperature hysteresis. 

Power is supplied to the controller from 120 Vac household power. The 120 Vac is 
converted to a lower value dc voltage and is supplied to the collector temperature detection 


circuitry and collector sensor. The collector temperature detection circuitry receives the 
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Figure 14. Block diagram of a basic differential temperature controller. 


collector temperature as an electrical resistance that changes with temperature. As a result 
of this resistance and de power supplied to the circuit, a voltage is developed within the 
detection circuit that is fed to the comparing circuitry. The dc power from the power 
supply is also supplied to the load temperature detection circuitry, which operates 
identically to the collector temperature detection circuitry. The two voltages fed to the 


comparing Circuitry are compared constantly. When the comparing circuitry sees that the 
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Figure 15. Differential temperature controller output wiring. 


collector temperature is greater than that of the load by a set amount, it senc . a signal to the 
output circuitry. ‘se circuitry will energize a relay (Figure 15), whose constants can be 
sent with whatever power is required to actuate pumps, valves, etc. When the output of the 
DTC is closure of a set of relay constants, the wiring of the output is as shown in Figure 


15. 


E. PHOTOVOLTAIC ARRAY DESIGN 

The photovoltaic or solar cell converts incoming light into direct current. Figure 16 
illustrates a construction of a typical photovoltaic cell which is very similar to a silicon 
diode. When the front and back contacts are connected to an electrical load, current will 
flow as a result of the induced voltage and dc electric energy will be generated. The solar 
cells are assembled into modules which are the smallest units that the consumer will 
consider. The voltage of a module is determined by the number of solar cells connected in 


series and usually each cell contributes 0.5 Volts. 











FRONT 
CONTACT 
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SEMICONDUCTOR 
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ELECTRIC 
LOAD BOTTOM 
SEMICONDUCTOR 





Figure 16. Photovoltaic (solar) cell construction. 


The grid - connected photovoltaic system, which is considered as the electric power 
source for the pump P, and the control unit, consists of an array and power conditioner 
interconnected with ac loads. The power conditiong unit (PCU) controls the operation of 
the array, starting it up in the morning, shutting it off at night, and maintaining the proper 


voltage for maximum power extraction. It also converts the variable - voltage dc electrical 


By) 





TN ee ee 


output of the array to constant voltage ac electric energy for the load. During a normal 


operation, the power conditioner controls the array voltage to track automatically maximum 





power whithin certain limits. When array current exceeds the capability of the power 
conditioner, such as may occur on a clear day, the unit will increase the array voltage to 
limit the excess current. When the array power falls below a minimum level the PCU will 
stop automatically converting energy and shutdown. 

The array consists of several series strings of modules which are then paralleled to 
achieve the desired power level. The series strings are designed to provide a voltage near 
the center of the power condotioner's input voltage operating range. Since the array 
configuration depends on PCU properties, the design process begins by choosing the 
power conditioner and then designing the array. To achieve a desired output the relation 
between the rating of the power conditioner and the array size is determined by the 
efficiency of the power conditioner and the average solar insolation. Since the usual PCU 
efficiency is 74%, then a dc input of about 1.35 kW is needed for each 1kW,, and is given 
by an array of 12 m? (130 ft2) approximately. In order to provide ac electric power 8.5 kW 
to the pump P, and the control unit, it is required an array with output dc power 
8.5kW/0.74 = 11.5 kW. For a formal module the manifacturer's specifications under 


standard operating conditions are the following: 


a. Voltage 12 V 
b. Current 2.5A 
c. Power 30 W 


Then the required area for the array to provide the desired ac power is: 


S we eeaiseis agai? 
array ~ “* kW 
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Figure 17. Interconnection of the photovoltaic system with the load. 


Finally determining the series / parallel arrangement of modules taking into account the 


above specifications,we can calculate the number of modules. It is given by: 


5 
sey TONED =, NOM 30s 3 dea wiodiiies 
Module Power 30 


while the number of modules per source circuit (series connection) is given by: 


39 








Array Voltage _ 120 _ 19 modules 
Module Voltage 12 


Then, we determine the number of modules in the array as 390 with 39 modules in 
parallel connection producing a peak current 97.5 A and 10 in series providing the nominal 
voltage of 120 Vdc to the pump and the control unit (Figure 17). 

Given the specifications of this particular module, the average cost of a panel with a 
power of 110 Watts is $ 500, that is, the cost of the array is $ 52,000 including the 


installation. 
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IV. CONCLUSIONS 


Economic consideration of the designed solar heating system is very important to 
examine and analyze its economic merit and feasibility. The alternative must be compared 
with the already existing fuel heating system on an equivalent economic basis. 

Using records and information from the boiler plant, we calculate the average amount 
of steam used to heat the pool. From boiler data, 4400 ft? of gas produce 4925 lbs of 
steam, so 1 ft? corresponds to 1.12 Ibs of steam or 1084 Btu. Also each therm costs $0.90 
and the average amount of steam is 20,000 Ibs per day. To calculate the cost per day, we 


use the equation: 


lft gas 1084Btu Itherm = $0.90 
CrP Aine Laon xX 
L121lbs steam = ft’ gas_~—s- 100000 Btu_—s therm 





Cost = (Ibs of steam used) x 


which gives a daily cost of $175 or $31,500 for a season of 180 days. 

For the designed system the capital cost will be $156,000 + $52,000 + $3,000 = 
$211,000. Assuming a service life of 10 years with an interest i = 5 % and using the 
straight - line depreciation plus average interest method, we approximate the true equivalent 


annual cost. These two variables are obtained from the equations: 


Straight - line depreciation = ——— (4.1) 





Average interest = (P - L(+] ae -)+ Li (4.2) 
2 n 











where, 

P = Capital cost (P = $211,000) 

L = Estimated salvage price (In our case we assume it will be 50 % of the capital cost) 

i = Interest rate 

n = Estimated life time 

The annual cost is the depreciation plus the average interest. From eq. (4.1), we geta 
depreciation of $21,100 and an average interest $8,176. Then the annual cost is $29,276. 

However, with our new designed system we are achieving a benefit of roughly $2,200 
per year. In this evaluation, the maintenance cost and the personnel payments, which are 
much lower for a solar heating system, were not taken into consideration. Its efficiency is 


also dependent on weather conditions, but not on fuel price variations. 
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TABLE A-2 THERMOPILYSICAL PROPERIES OF MASONRY 
MATERIALS. 


TYPICAL PROPERTIES AT 300 K 


DENSITY, THERMAL SPECIFIC 
p CONDUCTIVITY, k HEAT, c, 
DESCRIPTION /COMPOSITION (kg /ar) (W/m - K) (J/kg - K) 
Building Boards 
Asbestos-cement board 1,920 0.58 _ 
Gypsum or plaster board 800 0.17 _ 
Plywood 545 0.12 1,215 
Sheathing, regular density 290 0.055 1,300 
Acoustic tile 290 0.058 1,340 
Hardboard, siding 640 0.094 1,170 
Hardboard, high density 1,010 0.15 1,380 
Particle board, low density $90 0.078 1,300 
Particle board, high density 1,000 0.170 1,300 
Woods 
Hardwoods (oak, maple) 720 0.16 1,255 
Softwoods (fir, pine) 510 0.12 1,380 
Masonry Materials 
Cement mortar 1,860 0.72 780 
Brick, common 1,920 0.72 835 
Brick, face 2,083 13 _ 
Clay tile, hollow 
1 cell deep, 10 cm thick — 0.52 — 
3 cells deep, 30 cm thick — 0.69 _ 
Concrete block, 3 oval cores 
sand /gravel, 20 cm thick _ 1.0 _ 
cinder aggregate, 20 cm thick - 0.67 -- 
Concrete block, rectangular core 
2 cores, 20 cm thick, 16 kg _ 1.1 — 
same with filled cores - 0.60 _ 
Plastering Materials 
Cement plaster, sand aggregate 1,860 0.72 _ 
Gypsum plaster, sand aggregate 1,680 0.22 1,085 
Gypsum plaster, vermiculite 720 0.25 _ 


aggregate 





TABLE A-3 THERMOPITYSICAL PROPERIES OF THE AIR. 


T p c p- 10? »- 10° k- 10° a - 10° 


P 
(K) — (kg/m’) (kJ/kg K) (N-s/m?) — (m’/s) (W/m-K) — (m?/s) Pr 
Air 
100 3.5562 1.032 712 2.00 9.34 2.54 0.786 
150 2.3364 1.017 103.4 4.426 13.8 5.84 0.758 
200 1.7458 1.007 132.5 7.590 18.1 10.3 0.737 
250 1.3947 1.006 159.6 11.44 22.3 15.9 0.720 
300 1.1614 1.007 1846 15.89 26.3 22.5 0.707 
350 0.9950 1.009 208.2 20.92 30.0 29.9 0.700 
400 0.8711 1.014 230.1 26.41 33.8 38.3 0.690 
450 0.7740 1.021 250.7 32.39 37.3 47.2 0.686 
500 0.6964 1.030 270.1 38.79 40.7 56.7 0.684 
$50 0.6329 1.040 288.4 45.57 43.9 66.7 0.683 
600 0.5804 1.051 305.8 $2.69 46.9 76.9 0.685 
650 0.5356 1.063 322.5 60.21 49.7 87.3 0.690 
700 0.4975 1.075 338.8 68.10 §2.4 98.0 0.695 
750 0.4643 1.087 354.6 16.37 549 109 0.702 
800 0.4354 1.099 369.8 84.93 57.3 120 0.709 
850 0.4097 1.110 384.3 93.80 59.6 131 0.716 
900 0.3868 1.121 398.1 102.9 62.0 143 0.720 
950 0.3666 1.131 4113 112.2 64.3 155 0.723 
1000 0.3482 1.141 424.4 121.9 66.7 168 0.726 
1100 0.3166 1.159 449.0 141.8 71.5 195 0.728 
1200 0.2902 1.175 473.0 162.9 76.3 224 0.728 
1300 0.2679 1.189 496.0 185.1 82 238 0.719 
1400 0.2488 1.207 530 213 91 303 0.703 
1500 0.2322 1.230 557 240 100 350 0.685 
1600 0.2177 1.248 584 268 106 390 0.688 
1700 0.2049 1.267 611 298 113 435 0.655 
1800 0.1935 1.286 637 329 120 482 0.683 
1900 0.1833 1.307 663 362 128 $34 0.677 
2000 0.1741 1.337 689 396 137 589 0.672 
2100 0.1658 1.372 715 431 147 646 0.667 
2200 0.1582 1.417 740 468 160 714 0.655 
2300 0.1513 1.478 766 506 175 783 0.647 
2400 0.1448 1.558 792 547 196 869 0.630 
2500 0.1389 1.665 818 589 222 960 0.613 
3000 0.1135 2.726 955 841 486 1570 0.536 
Ammonia (NH) 
300 0.6894 2.158 101.5 147 24.7 16.6 0.887 
320 0.6448 2.170 109 16.9 27.2 19.4 0.870 
340 0.6059 2.192 116.5 19.2 29.3 22.1 0.872 
360 0.5716 2.221 124 21.7 31.6 24.9 0.872 
380 0.5410 2.254 131 242 340 27.9 0 869 
400 0.5136 2.287 138 26.9 37.0 31.5 0.853 
420 0.4888 2.322 145 29.7 40.4 356 0.833 
440 0.4664 2.357 152.5 32.7 43.5 39.6 0.826 
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TABLE A-5 EMISSIVITY OF NONMETALLIC MATERIALS. 








TEMPERATURE EMISSIVITY 





DESCRIPTION / COMPOSITION (K) € 
Silicon carbide (1) 600 0.87 
1000 0.87 
1500 0.85 
Skin (h) 300 0.95 
Snow (A) 273 0.82~0.90 
Soil (4) 300 0.93~0.96 
Rocks (h) 300 0.88-0.95 
Teflon {h) 300 0.85 
400 0.87 
500 0.92 
Vegetation (A) 300 G.92-0.96 
Water (h) 300 0.96 


“Adapted from Reference 1. 
"Adapted from References 1, 9, 20, and 21. 
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TABLE A-6 THERMAL AND VAPOR DIFFUSIVITY FOR DRY AND 
SATURATED MOIST AIR. 








Temperature Degree of D 
F Saturation Sfhr 








0 
0.901 
j 
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TABLE A-7 HEAD LOSS AND WATER VELOCITIES INTO PEPES. 


Velocity measured in {t./sec., Loss in feet of water head per 100 ft. of pipe. 


a yy," yy" 1” rye iy” 2" 2” 3” 3” 4” 
PE 
MUN. Vel Loss Vel Less Vel Low Wel Los Vel Less Wel Lose Wel Less Wel Less Wel Loss bel Loss 


2 274 OTF ban 199 
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APPENDIX B (SWIMMING POOL PLANS). 





Figure 18. Three-dimension plan of the Naval Postgraduate School 


swimming pool. 
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Figure 19. Heating system of the swimming pool. 
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APPENDIX C (OPTIMIZATION). 


echo on 





cle 

% {t plots the total solar irradiation (W/sq.m) vs tilting angle of the 

% collector for the winter time. 

x=[50,51,52,53,54,55,56,57,58,59,60,61,62]; 
y=[41685,41760,41823,41873,41910,41934,41946,41944,41930,41904,41864,.41812,4 
1747); 

plot(x,y); grid: 

xlabel(‘Tilting angle of the collector (degrees)'); 


ylabel('Total solar irradiation (W/sq.m)'); 


echo on 

cle 

% It plots the total solar irradiation (W/sq.m) vs tlting angle of the 

% collector for the summer time. 

x=[4,5,6,7,8,9,10,11,12,13,14,15,16]; 
y=[44234,44313,44378, 44431 ,44469,44494,44506,44504,44488 44459 444 16,44360,4 
4290]; 

plot(x,y); grid; 


xlabel(‘Tilting angle of the collector (degrees)’), 


ylabel(‘Total solar irradiation (W/sq.m)’); 








echo on 


cle 
% It plots the total solar irradiation (W/sq.m) vs tilting angle of the % %collector for all 
over the year. 
x=[27,28,29,30,31,32,33,34,35,36,37,38,39]; 
y=[78680,78806,78907,78985,79038,79067,79072,79053,79010,78943,78852,78737,7 
8598); 
plot(x,y); grid; 
ylabel(‘Total solar irradiation (W/sq.m)’); 
xlabel(‘Tilting angle of the collector (degrees)’); 
echo on 
cle 
a=0.0698 1317007; 
theta=0.9322967; 
Snav=0; 
Sntotal=0; 
$=725; 
for x=80:264; 
thetad=0.41015237422*sin(2*pi*(x-80)./365); 
phi=asin(-tan(thetad)./tan(theta)); 
Snav=(S/pi)*(cos(thetad)*sin(theta+a)*cos(phi)+sin(thetad) *cos(thetat+a)*((pi/2)- 
phi)); 
Sntotal=Sntotal+Snav; 
end 


Sntotal 
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